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1. Introduction 

The region surrounding the Galactic centre (GC) is among the brightest and most complex in 
high-energy 7 -rays, with on-going massive star formation providing all types of known or suspected 
cosmic ray (CR) and 7 -ray sources. The GC also houses a 10 6 M 0 black hole (e.g., [1]) and the 
region is predicted to be the brightest source of 7 -rays associated with annihilation or decay of 
massive weakly-interacting particles (see the reviews by, e.g., [2, 3, 4]). Despite detection in the 
100 MeV to GeV range by the EGRET instrument on the Compton Gamma-Ray Observatory [5] 
and at higher energies by the H.E.S.S. Cherenkov array [ 6 , 7] the characterisation of the 7 -ray 
emission for < 100 GeV energies in the region surrounding the GC has remained elusive. 

The 7 -ray emission in the Galaxy is predominantly due to the interactions of CR particles with 
the interstellar gas and radiation fields. This interstellar emission is a fore-/background against 
which 7 -ray point sources are detected. In the Galactic plane, and particularly toward the GC, the 
intensity of this emission makes disentangling the contributions by 7 -ray point sources and truly 
diffuse processes particularly challenging. 

Since 2008 the Large Area Telescope instrument on the Fermi Gamma-Ray Space Telescope 
(. Fermi-LAI ) has been taking data in the range 20 MeV to more than 300 GeV energies. In this 
contribution, an analysis is described of the 7 -ray emission observed by the Fermi-LAI during the 
first 62 months of the mission toward the inner Milky Way that characterises the 15° x 15° region 
in Galactic coordinates centred on the GC. This encompasses the innermost ~ 1 kpc where the CR 
intensities, interstellar gas and radiation field densities are highest but most uncertain, and signa¬ 
tures of new physics may be detectable. The analysis uses multiple interstellar emission models 
(IEMs) together with an iterative fitting procedure to self-consistently determine the contributions 
by diffuse and discrete sources of high-energy 7 -ray emission. The GALPROP CR propagation 
code 1 (e.g., [ 8 , 9]) is used to calculate components of IEMs that are fit to the Fermi-LAI data to 
predict the interstellar emission fore-/background toward the 15° x 15° region. Candidate locations 
of point sources are found using a wavelet-based algorithm [10, 11]. These are used together with 
the IEMs to define a model for the emission of the region, which is then optimised in a maximum- 
likelihood fit to determine the contribution by CR-induced diffuse emission from the innermost 
~ 1 kpc and 7 -ray point sources. 

2. Methodology 

2.1 Data Selection and Preparation 

The analysis employs events with reconstructed energy in the range 1-100 GeV, where the 
effective area of the LAT is largest and not strongly dependent on energy. To allow the best separa¬ 
tion between point sources and the structured interstellar emission in the analysis procedure, only 
front-converting events are used. 

Events are selected from approximately 62 months of data from 2008-08-11 until 2013-10-15 
using the standard low-residual CR background “Clean” events from the Pass 7 event selections 2 . 

^or a detailed description of the GALPROP code the reader is referred to the dedicated website: 
http ://galprop. stanford.edu 

2 The reprocessed data and instrument response functions P7REP_CLEAN_V15 are employed. 
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Exposure maps and the PSF for the pointing history of the observations were generated using the 
standard Fermi-LAI ScienceTools package (version 09-34-02) available from the Fermi Science 
Support Center 3 . 

2.2 Fore-/Background Modelling 

Specialised IEMs are constructed to estimate the fore-/background interstellar emission to¬ 
ward and through the 15° x 15° region about the GC. The results of the Fermi-LAI team study 
[12], which compared GALPROP- generated IEMs normalised to local CR measurements with y- 
ray data, are used. A major uncertainty affecting predictions of the interstellar emission toward the 
inner Galaxy is the spatial distribution of CR sources. The pulsar distribution (“Pulsars” [15]) and 
the distribution of OB stars (“OB stars” [16]) employed in the aforementioned Fermi-LAI study 
encapsulate this because they represent reasonable extremes for the Galactocentric radial depen¬ 
dence 4 . These CR source distributions are used 5 to calculate “baseline” y-ray intensity templates 
for Galactocentric radial annuli for the standard processes (7T°-decay, IC, Bremsstrahlung 6 ). The 
annular intensity maps are used as templates together with an isotropic component and a model for 
y-ray emission associated with the Loop I supernova remnant (SNR) employing a two-component 
spatial template from [13] with a power-law spectral model for each, and point sources from the 
Third Fermi source catalogue (3FGL; [14]). The combined model for each IEM is fit to the Fermi- 
LAI data excluding the 15° x 15° region about the GC using a maximum-likelihood method. 

Figure 1 shows the fractional residuals, ( data — model)/model, for the Pulsars IEM for the 
1-3.16 GeV energy band for the (left to right) baseline, and two IEM variants of this model that 
have been scaled to the data (described below). The regions not used in the fitting procedure are 
explicitly masked in the figure. They are not used because of localised extended excesses that are 
most likely unrelated to the large-scale interstellar emission. 

Outside of the Galactic plane the fractional residuals are substantially reduced for the IEM 
shown in the central panel compared to the baseline. This IEM variant is termed ‘intensity-scaled’ 
because only the normalisations of the individual intensity maps are adjusted in the scaling proce¬ 
dure. The lower residuals for this IEM are due to the scaling of the 7T°-decay interstellar emission 
for the local annulus, and of the IC component generally. Along the Galactic plane the > 30% 
under-prediction by the baseline model is reduced to < ±10% after this scaling, except for scattered 
regions. However, at higher energies (not shown) the intensity-scaled IEM tends to over-predict 
the data along the plane interior to the solar circle. To account for this, another variant is developed 
that includes additional degrees of freedom to the spectrum of the 7T°-decay emission for annuli 
iterior to the solar circle. It is fit to the data following the same procedure as the intensity-scaled 
IEM. This variant is termed ‘index-scaled’ and has a fractional residual < ±10% for longitudes 
Z ^ — (15 — 70)° for the 1-3.16 GeV band (Fig. 1, right panel), with a slight increase in the residual 

3 http ://fermi .gsfc.nasa. go v/s sc/data/analy si s/ 

4 The reader is referred to Fig. 1 in [12] for a visualisation of the Galactocentric radial variation. 

5 The models assume an axisymmetric cylindrical geometry for the CR confinement volume with a halo height 
Zh = 6 kpc and maximum radial boundary R[ 2 = 30 kpc, and all other parameters the same so that the only difference is 
the assumed spatial distribution for the CR sources. 

6 Bremsstrahlung is a minor component over the 1-100 GeV energy range used in the analysis, and is held constant 
at the GALPROP predictions for each IEM. 
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Figure 1: Residual fractional counts [data —model)/model in the 1-3.16 GeV energy range for the baseline 
Pulsars model (left), intensity-scaled Pulsars model (centre), and index-scaled Pulsars model (right) fitted 
following the procedure described in the text. The baseline model does not include a model for the Loop-I 
SNR, resulting in the positive residuals in the northern Galactic hemisphere inside the solar circle. The 
maps are calculated for a HEALPix [17] order 8 pixelisation (~ 0.23° resolution) and smoothed with a 1° 
FWHM Gaussian. The positive residual at mid-to-high latitudes interior to the solar circle is due to mismatch 
between the data and the relatively simple Loop I model. The residuals close to the plane caused by this 
simple model are lower and do not affect the analysis of the 15° x 15° region about the GC. 

for the corresponding positive longitude range; the residuals improve to a similar degree at higher 
energies. Qualitatively, similar results for the scaled OB stars IEMs (intensity/index-scaled) are 
also obtained. It is not straightforward to identify a best IEM after fitting because the qualitative 
improvement for each over the corresponding baseline IEM is similar. Consequently all 4 IEMs 
are used to estimate the fore-/background below. 

2.3 Modelling 15° x 15° Region about the Galactic Centre 

Modelling of the interstellar emission and point sources over the 15° x 15° region about the 
GC is accomplished using an iterative procedure that identifies point-candidates (‘seeds’) and fits 
for their fluxes and spectra together with the interstellar emission from 7T°-decay associated with 
the neutral gas and IC components for the innermost annulus using a maximum-likelihood method, 
while the fore-/background interstellar emission (determined above) outside the inner ~ 1 kpc is 
held constant. 

Point-source seeds are identified using the wavelet analysis algorithm PGWave [10, 11]. Ap¬ 
plication of PGWave identifies true point sources, as well as structures in the interstellar emission 
that are indistinguishable from point sources due to the finite angular resolution and statistics of 
the Fermi-LAI data, without dependence on the specifics of an IEM. The spectra of candidates are 
initially evaluated using PointLike , a package for maximum-likelihood analysis of Fermi-LAI data 
[18, 19]. This package is also used to optimise the positions of seeds from the PGWave- determined 
list. For the only extended source that has been previously identified in the region, the W28 su¬ 
pernova remnant [20], the spatial template and spectral model employed for the 3FGL analysis are 
used. Its spectral parameters are refit during the maximum-likelihood procedure. 

Point-source candidates are combined with their PointLike trial spectra together with the fore- 
/background models in a second maximum-likelihood fit. A binned likelihood fit is performed 
using the Fermi ScienceTool gtlike. The templates for the 7T°-decay related 7 -ray intensity from 
H I and CO, and the IC emission, in annulus 1 are freely scaled in the fitting procedure. The 
point-source seed detection is run again on the residual maps to find fainter sources that were 
missed in the first iteration. PointLike is again used to determine their initial spectra and optimise 
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their localisations. The combined point-source candidate list from the first and second iterations 
7 , and the interstellar emission components for annulus 1 are fit using gtlike. The results of the 
maximum-likelihood fit are values and confidence ranges for the coefficients of the H I annulus 1, 
CO annulus 1, IC annulus 1, as well as the TS , fluxes and spectra for the point sources. All point 
sources with a maximum-likelihood determined TS > 9 are included in the model. 


3. Results 

3.1 Interstellar Emission 


Preliminary Preliminary 




Figure 2: Differential flux components obtained for the Pulsars intensity-scaled (left) and index-scaled 
(right) IEMs for the 15° x 15° region about the GC. Line styles, as in figure legend. Solid circles: data. 

Figure 2 shows the differential spectra of the individual components obtained for the Pulsars 
intensity-scaled and index-scaled IEMs integrated over the 15° x 15° region about the GC (the 
results for the OB stars IEMs are similar, but not shown due to space limitations). The figure 
separates the emission components in terms of the contributions by 7T°-decay and IC for annulus 1, 
the interstellar emission fore-/background, and point sources over the region. As expected, the fore- 
/background dominates for each IEM, which is predominantly 7T°-decay in origin. Interestingly, the 
IC is the dominant interstellar emission component over the inner ^1.5 kpc, and is much higher 
than predicted by GALPROP for the baseline IEMs. Its contribution to the total flux depends on 
the IEM and whether the residual is fitted (Section 3.3). The variation of the IC flux for annulus 1 
over all 4 IEMs is within a factor ^ 1.5. The interstellar emission from the neutral gas 7T°-decay 
is considerably suppressed compared to the GALPROP predictions. For both Pulsars and OB stars 
models the interstellar emission fore-/background is harder at high energies for the index-scaled 
IEMs, which reduces the spectral residuals for energies >10 GeV. 
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Figure 3: Point sources for 3FGL (left panel) and 1FIG (right panel, for Pulsars intensity-scaled IEM) 
overlaid on the total counts for the 15° x 15° region about the GC. Left panel symbol key: filled squares, 
‘flagged’ 3FGL sources; filled triangles, other 3FGL sources; upright crosses, 3FGL sources with a multi¬ 
wavelength association. Right panel symbol key: filled circles, 1FIG sources with TS > 25; angled crosses, 
1FIG source candidates with TS < 25; upright crosses, as in left panel. Colour scale is in counts per 0.05 2 
degree pixel. 


3.2 Point Sources 

The analysis finds 48 point sources over the 15° x 15° region about the GC with a TS > 25 
for the Pulsars intensity-scaled IEM, which was the model used for the point-source positions and 
localisation uncertainties 7 8 . This collection is denoted the 1 st Fermi Inner Galaxy (1FIG) source list. 
Figure 3 shows the point sources from the 3FGL and 1FIG overlaid on the total photon counts for 
the 15° x 15° region about the GC. The 3FGL sources are separated according to whether they have 
an analysis flag set in the 3FGL catalogue: flagged sources indicate their properties depend on the 
IEM or other details of the analysis in the region. The density of flagged 3FGL sources is higher 
out of the Galactic plane than that of the 1FIG sources, even if the TS < 25 source candidates are 
included. This can be partly attributed to differences in the IEMs used for the respective analyses. 

Some of the 1FIG sources and source candidates are likely due to mismodelling of the inter¬ 
stellar emission. However, it is difficult to determine the precise fraction. Multi-wavelength asso¬ 
ciation with objects of known y-ray emitting source classes can suggest counterparts. There are 14 
1FIG source and source candidate overlaps when comparing with SNRs from Green’s SNR cata¬ 
logue 9 [21] and pulsars from the ATNF catalogue 10 [22]. Multiple overlaps occur across and within 
the catalogues, e.g., SNR 354.1+00.1 and PSR J1701-3006A,D,E overlap with 1FIG J1701.1-3004. 
Three of the 1FIG sources overlap with SNRs that have not previously been detected in high-energy 
y-rays; follow-on studies are required to better characterise their spatial and spectral properties. 

7 No new significant excesses that could be identified as viable seeds were found after two iterations of the procedure. 

8 The optimisation of the source seed localisation does not have a strong dependence on the IEM. 

9 http://w w w. mr ao. cam. ac. uk/ survey s/snr s / 

10 http://www.atnf.csiro.au/people/pulsar/psrcat/ 
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3.3 Residuals 

Preliminary 



5 0-5 


Preliminary 



Preliminary 



5 0-5 


Longitude (degrees) 


Longitude (degrees) 


Longitude (degrees) 


Figure 4: Longitude profiles for 1-1.6 (left), 1.6-10 (middle), and 10-100 GeV energies (left), re¬ 
spectively, of the residual counts, data — model , for the Pulsars index-scaled IEM after fitting for in¬ 
terstellar emission and point sources across the 15° x 15° region. Line styles: black/solid, total model; 
cyan/dashed, fore/background interstellar emission; green/dotted, point sources; magenta/solid, IC from an¬ 
nulus 1; blue/solid, 7T°-decay from annulus 1. Point styles: red, data; black, residual counts. The lower 
sub-panel for each profile gives the residual counts after the model has been subtracted from the data. The 
error bars are statistical. Profiles for the residuals counts for other IEMs display similar features with the 
major difference being the number of counts. 

Figure 4 shows the longitude profiles for the 1-1.6, 1.6-10, and > 10 GeV energy ranges 11 for 
the Pulsars index-scaled model, which is the model that has the lowest fractional residual across 
the 1-100 GeV energy range . The lower sub-panel for each figure gives the residual counts, 
data — model. (The features are mostly the same for the other IEMs.) While there is considerable 
statistical noise, there is some indication that the residual counts are distributed asymmetrically in 
longitude about the GC below 10 GeV. 

The model over prediction at the lowest energies is primarily correlated with the Galactic 
plane, which could be due to mismodelling of the gas component of the IEMs. Some of the posi¬ 
tive residual in the few GeV range could be due to an extended component that is more concentrated 
toward the GC compared to the IEM components. However, it is difficult to establish properties 
for an additional component not presently included in the model for the region. Different spatial 
models are considered for the positive residuals: a squared Navarro-Frenk-White (NFW) profile 
[23] for dark matter annihilation, 2D Gaussians with HWHM 1°, 2°, 5°, and 10°, and the CO annu¬ 
lus 1 gas template smoothed with a 2° Gaussian as a tracer of unresolved pulsars. An exponential 
cutoff power law is used for the spectral model, which has some flexibility to model a pulsar or a 
DM annihilation spectrum without supposing specific scenarios. The positive residual model is fit 
together with the interstellar emission and point sources (Sec. 2.3). Its flux and spectral parameters 
depend strongly on the interstellar emission fore-/background, and are also covariant with the nor¬ 
malisations of the standard interstellar emission processes (;r 0 -decay and IC) for annulus 1. Note 

11 Each band covers the energy intervals where the under/over-predictions for the fractional residuals are more promi¬ 
nent in Fig. 2. 


7 
































Observations of High-Energy y-rays Toward the Galactic Centre with the Fermi-LAT 


Troy A. Porter 


that not all of the positive residual is accounted for by the centrally peaked models that are consid¬ 
ered here. The remainder is distributed about the Galactic plane, but with no obvious correlation to 
the spatial templates used for the analysis. 
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